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We present a new approach to treat the perturbative QCD evolution of TMDs by resummation. 
We obtain a universal evolution kernel for the eight quark distributions while large logarithms are 
resummed up to next-to-next-to leading logarithms (NNLL). Contrary to earlier works this resum- 
mation gives unambiguous predictions avoiding the need to introduce an ad-hoc cut-off parameters 
and thus highly reducing the model-dependence. Applications to single-spin asymmetries, as the 
Sivers function, and unpolarized transverse-momentum-dependent parton distribution functions are 
given. The results agree very well with phenomenology. 
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Introduction: To formulate a complete, consistent 
theory of TMDs is one of the great challenges in modern 
QCD. TMDs are needed for all processes for which intrin- 
sic transverse parton momenta are relevant, which form 
a large group. For example spin-dependent transverse 
momentum asymmetries provide unique clues to clarify 
the internal spin, angular momentum and 3-dimensional 
structure of hadrons. While operating at different ener- 
gies, experiments and facilities such as HERMES, COM- 
PASS, JLab and BNL are pursuing intensive research 
programs to explore TMDs. On the theoretical side, 
Sivers [1] and Collins [2] asymmetries have been intensely 
studied (see [3] for a review of TMD functions in spin- 
physics), and have attracted much attention recently [4- 
10] in context with fundamental theoretical concepts. 
Basically, some of the observed spin-asymmetries are 
linked to the presence of gauge links in non-collinear non- 
local correlators needed to maintain gauge invariance. 
A complete formalism for the treatment of TMDs is 
still being developed. In this paper we refer to two ap- 
proaches, namely the one developed by three of us in [11] 
and the one developed mainly by Collins [12]. For certain 
choices of parameter(s) one can argue that the two are 
equivalent, despite the technical differences. Here we ad- 
dress a long-standing problem of mo del- dependence and 
show how it is resolved by a new resummation technique 
which relies, in part, on the extracted anomalous dimen- 
sion of the unpolarized quark-TMD up to 0{al), as was 
explained in [11]. 

The study of the unpolarized TMDs was pioneered by 
Collins and Soper [13, 14]. Collins' approach to TMDs 
consists of properly defining those quantities in a way 



consistent with a generic factorization theorem, extract- 
ing their anomalous dimensions and their evolution prop- 
erties. This approach relies mainly on taking some of the 
Wilson lines in the soft factor off-the-light-cone. When 
doing so, one introduces an auxiliary parameter C spec- 
ifying the measure of "off-the-light-coness" . A differen- 
tial equation with respect to Cj the Collins-Soper evo- 
lution equation is then derived and solved which resums 
large logarithms and determines the evolution of the non- 
perturbative TMDs with energy. 

Though Collins' approach is theoretically robust, two 
main drawbacks of it affect its phenomenological appli- 
cations. First is the difficulty to perform off-the-light- 
cone calculations (beyond 0{as)) compared to the on- 
the-light-cone ones. This, in practice, limits the accu- 
racy to which large logarithms can be resummed. Sec- 
ond, and in order to avoid the Landau pole, an ad hoc 
non-perturbative model has to be introduced, which de- 
pends on an arbitrary parameter fomax that can only be 
extracted through fitting to experimental data. It turns 
out that the choice of 6niax affects the evolution of TMDs 
not only at very low energies but at intermediate ones 
as well. This can be clearly noticed from the plots intro- 
duced below. 

In this letter we present a new resummation technique 
based on the formalism developed in [11]. With what we 
call the "resummed D" , explained below, we avoid the 
issue of Landau pole in a way where there is no need to 
introduce any ad hoc (cut-off) parameter, thereby elimi- 
nating the model-dependence to a precision far exceeding 
all current practical purposes. The parameter-free pre- 
diction of our formalism agrees perfectly with the results 



2 



obtained for the phenomenologically preferred value of 
bmax [15]. We see this as an important step towards a 
complete theory of TMDs. 

Definition of Quark-TMDs: Extending the work 
done in [11], we define a quark-TMD of a polarized 
hadron, coUinear with the +z direction with momentum 
P and spin S as 



(27r)3 



X $^^(0+,r-,fi)v/5(0+,0-,rl). (1) 



^ap stands for a purely coUinear matrix element, i.e., a 
matrix element which has no overlap with the soft region 
[16], and it is given by the bilocal correlator 

^Ip = {PS\ [CnaW^J] (0+,y-,y±) [WP^^p] (0) \PS} 



(2) 



relation between the anomalous dimensions of F and H, 
i.e., jF = --^IH, where = 2rcusp ln(Q2//i2) + 2jv 
is known at 3-loops level [21, 22]. This crucial result 
can be automatically extended to all TMDs defined in 
Eq. (1), since the anomalous dimension is independent 
of spin structure. 

Evolution Kernel: For spin-dependent TMDs the 
operator product expansion for coUinear PDFs fails. One 
has to resort to non-perturbative models, fitted at low 
energies to experimental data, to make predictions for 
higher energy experimental probes. Obviously, knowing 
the evolution kernel to as high accuracy as possible is 
key. 

Starting from Eq. (1) the evolution of a generic quark- 
TMD from an initial scale Qi to a final scale Qf is given 
by F{x,b;Qf) = F{x,b;Qi) R{b;Qi,Q f) where the evo- 
lution kernel .R is [8, 11, 12] 



The soft function S which encodes soft-gluon emission is 
given by 



S = 



iTr 
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(0)|0). (3) 



We should mention that Fap is free from rapidity diver- 
gences, thus it is independent of any parameter like C^. 

One can obtain the eight leading-twist quark-TMDs 
[17, 18], represented generically by F below, simply by 
tracing Fap with the Dirac structures ^, and "^^^ 
for unpolarized, linearly polarized and transversely polar- 
ized quarks inside a polarized hadron. The superscript T 
indicates transverse gauge-links Tn(n)-: necessary to ren- 
der the matrix elements gauge- invariant [19, 20]. The 
definitions of coUinear (Wn(n))i soft {Sn(fi)) and trans- 
verse {Tn(n)) Wilson lines for DY and DIS kinematics 
can be found in [11]. 

In [11] the anomalous dimension of the unpolarized 
TMD was given up to 3-loop order based on a factor- 
ization theorem for g^-dependent observables in a Drell- 
Yan process. In impact parameter space such factoriza- 
tion theorem for the hadronic tensor can be schematically 
written as: W = HFf/pFj^p, where H is the hard co- 
efficient which encodes the physics at the probing scale 
Q. This quantity is built, order by order in perturbation 
theory, by considering virtual Feynman diagrams only, 
i.e., no real gluon emission has to be considered. More- 
over, the quantity H has to be free from infrared physics, 
no matter how the latter is regularized. This should be 
the case whether one works on or off-the-light-cone. 

Since the factorization theorem given above holds, at 
leading twist, also for spin-dependent observables, one 
can apply the same arguments, based on renormalization 
group invariance, as for the unpolarized case, to get a 
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The D term can be obtained by noticing that the renor- 
malized F has to be well-defined when its partonic ver- 
sion is calculated pertubatively. This means that all di- 
vergences, other than the genuine long-distance ones has 
to cancel. This fundamental statement - that rapidity 
divergences cancel when the coUinear and soft matrix el- 
ements are combined according to Eq. (1)- allows one 
to extract all the Q^-dependence from the TMDs and 
exponentiate it with the function D [11] thereby sum- 
ming large logarithms of the form ln{Q^/q^). Apply- 
ing the renormalization group invariance to W we get 
dD/dlnfi = Fcusp by which we extract, for all the quark 
spin-dependent TMDs the function D at NNLO from the 
known cusp anomalous dimension at three- loops [22]. In 
this way we do not have to solve the Collins-Soper equa- 
tion. 

The function D contains large logarithms of the form 

2.2 w ^ 

L± = In 4^-2TE that have to be resummed. Matching 
the expansion D{b,fi) = X^^i dn{L_L) (fj)" to the ex- 
pansions rcusp(as) Tn-i (^)" and /3{as) = 
—'2asJ2^=il^n-i (j^) one gets the following recursive 
differential equation 



d'niLj 



ml3n-l-mdm{L±) , (5) 



where — ddn/dL±. Solving this equation we can ex- 
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press the coefficients d„ in tfie following way, 



Within Collins approach the evolution kernel is [5, 12] 



2d„ = ££iMi):+(^„L^)"-i 

Po n 
Po / 



^"-2 ' (n-l)rf0|„>2 + (n-l)-f|„>3 

^Po 



— Sri|n>4 + ^3-rn|ri>5 + -TTJ^i.^ ~ 6)\n>i 
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where 
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(l-n)i7f_^i+n+l 



1 

" 2 

-(n - l)(V'(7i) + 7B " 2)(?/'(n) + 7b)] 



(7) 



(r) 

Hn is the r-th order Harmonic Number function of n 
and "(/'(n) is the digamma function of n. We can perform 
the resummation of large logs in the D when a^Lj^ ^ 1 
and obtain with a = as/(47r) and X = a^oL^ 



2D = --^Hl-X) 
Po 



AFo X + ln(l - X) 
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For the last term in this expression we have found an an- 
alytic form of the sum of the series by using the approxi- 



mation H, 



(1) 



-1 -liiH+7£;-r2;i-T2;;^-rT2o;jT 

which is precise enough for our purposes. 

This resummed series is valid for \X\ < 1, however 
when analytically continued through Borel-summation 
its validity is extended to X — >■ — oo, which corre- 
sponds to 6 —7- 0. The maximum value of b where 
this series can be used, corresponding to X = 1, is 
be = 2e-T^//iexp[27r/(/3oas(Ai))], and takes 8.42 GeV"^ 
for M = GeV. 

For all practical applications, the resummed D will de- 
termine the evolution kernel, R, for all polarized and un- 
polarized TMDs. Moreover we are able to resum large 
logarithms up to NNLL accuracy. 

Next, we briefly present Collins approach and illus- 
trate the differences with ours by considering unpolarized 
quark TMD and Sivers function. 



2n ' 



1 

T2n^ 



1 



xexp-^ / —lF[o^s{^i)M^^ 




A'(b*,x/cr) 



(9) 



where the Collins-Soper kernel K is resummed through 
its RG-equation 



K 



Kib*;tib)- -^IK- (10) 



Large logarithms are cancelled by choosing /i^ = 
2g-7-B/5* [23]. To avoid the Landau pole that appears 
when Fourier transforming back to momentum space, one 
writes b* ^ b/ y^l + (6/5niax)^ instead of b. 6max is an ar- 
bitrary cutoff that accounts for the separation between 
perturbative and non-perturbative regimes. Below we 
have implemented the BLNY model to compare our ap- 
proach with Collins' with the following: gK{b) = ^b^ , 
for &„iax = 0.5 GeV"^ [24] and 



and 52 = 0.68 GeV^ 

t2 



52 = 0.184 GeV^ for 6„,ax = 1-5 GeV~' [15]. 

Within Collins' approach one usually identifies Q = 
and Cf ~ Q^f- With this choice one gets at 0{as): ~ 
jp and K ~ —2D, and those relations should hold to all 
orders in perturbation theory. However, 7^ beyond one- 
loop is not known and its determination would require 
an extremely difhcult perturbative calculations. 

Results and Conclusions: In Fig. 1 we show the 
evolution kernel from Qi = ^2.4 GeV to Q/ = 5 GeV 
with Collins' and our approaches. The impact of the 
differences on the evolution of two TMDs is illustrated in 
Fig. 2. 
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Evolution Kernel 

fo„,a,=0.5GeV-' 

fo„„=1.5GeV-' 

________ Resummed D at NLL 

Resummed D at NNLL 




Qi = ^J 2.4 GeV 
Qf = 5 GeV 



FIG. 1: Evolution kernel from Qi — v2.4GeV up to Qf 
5GeV and Qf — 91.19GeV with different approaches. 
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Evolved unpolarized TMDPDF 



i„„=0.5GeV' 

Resummed D at NLL 
Resummed D at NNLL 
Initial model 




Qi = y 2.4 GeV 
Qf = 5 GeV 
.v = 0.1 

U.J ] .U [ .J ^.0 ^ 

Evolved Sivers function (Bochum) 

t„„=0.5GeV"' 
i„„=1.5GeV"' 
Resummed D at NLL 
Resummed D at NNLL 
Initial model 




2.4 GeV 
= 5 GeV 
x = 0.1 



Evolved Sivers function (Torino) 




i,„„=0.5GeV-' 
fc,„„=I.5GeV"' 
Resummed D at NLL 
Resummed D at NNLL 
Initial model 



Qr 



2.4 GeV 
= 5 GeV 
.1- = 0. 1 



FIG. 2: Up quark unpolarized TMD and Sivers function 
(Bochum and Torino fits shown with opposite sign) evolved 
from Qi — V2.4GeV up to Qf = 5GeV with different ap- 
proaches to the evolution kernel. 



In order to perform the resummation of large loga- 
rithms consistently up to N*LL order (or N*~^LO in 
RG-improvcd perturbation theory) one needs the input 
shown in table I. In our approach one takes the resummed 
series in Eq. (8) up to the corresponding order i. In 
[4, 5, 8] 7/f was not implemented at 2-loop order, as it 
should be to get a complete NLL result. 



eled as a Gaussian: /jjp/p(a;; Qi) exp[— cr6y] with a = 
0.38/4 GeV^ [25], and we have taken the MSTW data 
set for the integrated PDFs [26]. The Sivers function 
at low energy is modeled following what are called the 
"Bochum" [27] and "Torino" [28] fits in [5]. 

The model-dependence within Collins' approach man- 
ifests itself as dependence on &niax and the fit parame- 
ters of, for example, the BLNY model. However within 
our approach this model-dependence is almost absent. In 
principle one would need a model from be onwards, where 
the resummed D is no more valid. But the value at which 
the mo del- dependence would enter has been shifted to 
larger values of b by the resummation, such that the re- 
maining model-dependence is negligible. Since the evolu- 
tion kernel and the TMDs that are evolved are very small 
in that region, one can set the kernel to zero with no ap- 
preciable effect on evolution. Thus no non-perturbative 
model is needed. 

As can be seen in Fig. 2, our implementation of Collins 
approach is consistent with [4, 5]. From this figure it is 
clear that if one uses Collins' approach rather than ours, 
the value of 6max has to be close to 1.5 GeV~^ as was 
indeed found in [15] by a comparison with experimental 
data. 

The definition of quark-TMDs and the new approach 
to determine the evolution kernel can be extended to 
gluon-TMDs [29] , which can be defined also on-the-light- 
cone. This new resummation technique can be applied 
as well to the evolution kernel of the complete hadronic 
tensor W (built with two TMDs). One major applica- 
tion would be to relate the ^T-spectrum of the Higgs 
boson production at the LHC at two different energies in 
a model-independent way. 
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Order 


Accuracy ~ a"L* 


7v Pcusp D 


N'LL 


n + l-i<k<2n {ai~^) 


... i „ * + l „ i 



TABLE I: Approximation schemes for the resummed TMD, 
where L — ln(Q^/^^) and al indicates the order of the per- 
turbative expansion. 

The unpolarized quark-TMD at low energy is mod- 
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